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Abstract. Much is known in the analysis of a finitely ramified self-similar fractal when 
the fractal has a harmonic structure: a Dirichlet form which respects the self-similarity of 
a fractal. What is still an open question is when such structure exists in general. In this 
paper, we introduce two fractals, the fractalina and the pillow, and compute their resistance 
scaling factor. This is the factor which dictates how the Dirichlet form scales with the self- 
similarity of the fractal. By knowing this factor one can compute the harmonic structure 
on the fractal. The fractalina has scaling factor (3 + \/4T)/16, and the pillow fractal has 
scaling factor 



1. Introduction 

The first step in understanding analysis and probability on a fractal space is calculating 
a Dirichlet form on a suitable domain of functions. As seen in IFOTl 11 . a Dirichlet form is 
equivalent to a Laplacian operator, and to a diffusion process on the fractal. The standard 
references for the theory of analysis and probability on fractals are fKigOT IStr06l IBN98L 
However, techniques used in computing Dirichlet forms are often algebraic — the current 
work rewrites the problem in terms of solving equations of homogeneous rational functions. 
This is the first step in using the theory of curves and surfaces to classify the space of forms 
on a given fractal. 

We calculate the Dirichlet forms on two self-similar fractals, the modified fractalina 
Sierpinski gasket and the pillow fractal. Recently there has been much progress in estab- 
lishing the existence and uniqueness of these forms, for example I BSU08, CS071IBCF"'"07[ 
IHMT06IIPei081 . Despite this, there are still simple examples, such as the fractalina, which 
elude understanding. In IINT08 1. the theory of self-similar groups is used to establish the 
existence and uniqueness of a scaling factor for the pillow fractal. The current work follows 
up by determining the exact value of this factor. 

Self-similarity — the quality of looking the same on all scales — is one of the most 
important concepts in the study fractal geometry and analysis on fractal spaces. A compact 
metric space K is called self-similar if there exists a finite collection of injective contraction 
maps {i^i}i=i with the property that 

N 
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Analysis on these fractals relies on Dirichlet forms. The classical example of a Dirichlet 
form is the energy integral / Vu • on W^. Dirichlet forms are bilinear forms the ex- 
istence which is equivalent to Laplacian operators and diffusion processes. For more on 
this relation, see IFOTllL A Dirichlet form on a self-similar fractal is self-similar if for 
some Pi > 0, i = 1, . . . , N 

N 

(1.1) <f(w,u) = ^p,^(uoF„woF,). 

i=l 

In the case where all of the pi = p are equal, we call p^^ a resistance scaling factor. Cal- 
culating this factor allows us to much better understand the energy form on the limiting 
fractal. We calculate this factor for two specific self-similar structures, the modified frac- 
talina Sierpinski gasket, or fractalina for short, and the pillow fractal. 

One particularly nice property of these fractals is that they are post-critically finite (p.c.f.) 
self-similar structures — there is a finite set Vq, called the fractal boundary, such that 
Fi{K) n Fj{K) C Fi{Vo) n -^,(^0), for all i j. A theory of analysis exists for a 
p.c.f. self-similar structure with a harmonic structure, see | [Kig 93 Ki gOT][Str06l . This liter- 
ature also details the existence of a self-similar Dirichlet form for highly symmetric nested 
fractals. 

The exact conditions for existence and uniqueness of a harmonic structure on self-similar 
fractals is still not fully understood, even for p.c.f. self-similar structures. There has been 
much progress in recent years, see liMet03l[Sab971 . In IIHMT06L the symmetry conditions 
are greatly relaxed. The existence of harmonic structures on other large classes of fractals 
is found in |Pei08|, including generalized Vicsek sets. 

The works |.CS07i IBSU08II identify a class of energy forms on non-equilateral Sierpin- 
ski gaskets. These energy forms are indexed by a 2 dimensional smooth manifold. This is 
proven by showing that the self-similar scaling factors of the forms satisfy certain polyno- 
mial relations and are thus contained in an algebraic variety. Even though the energy forms 
are not unique, algebraic techniques were used to classify them. 

Fractalina is a modified version of the Sierpinski gasket. Like the gasket, fractalina is 
generated by three contraction mappings, but one of the contraction features a 180° rotation. 
This rotation introduces a few complications. The boundary of this fractal now consists of 
4 points, and the fixed point of the rotated contraction is no longer in the boundary of the 
fractal. In section|5]we introduce the fractalina and compute its resistance scaling factor to 
be (3 + \/4T)/16. 

In section [6] we introduce the pillow and calculate the scaling factor ^^2. The solution 
is in interesting result because in most of the known examples, the resistance scaling 
factor is the solution to a quadratic equation. Establishing a scaling factor also proves that 
there is a self-similar resistance form on both of these fractals. 

In section[2]we introduce rigorous definitions of self-similar structures, and discuss how 
to approximate energy forms on these structures by looking at finite subsets. In section [3] 
we reduce the problem of defining a self-similar Dirichlet form to that of solving a matrix 
equation. Thus the problem becomes a matter of solving a finite number of rational equa- 
tions, a problem which is accessible to computer algebra systems. Section |4] discusses the 
techniques from electrical circuits which we use in the paper. 

Acknowledgements: The authors would like to thank Alexander Teplyaev for his insight 
and guidance. 
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2. SELF-SIMILAR STRUCTURES 



This section describes that basics of the theory of self-similar structures and their Dirich- 
let forms. For a more extensive treatment, see | Kig01| . 



Consider the set X = {!,... ,N}, which we shall call the alphabet. Then X" = 
{wiW2 ■ ■ ■ Wn \ Wi G X} will be the words of length n. We take X* — [J^^i X" to be 
the collection of all finite words, and S ^ YVuLi ^ '^o '^^e set of all right-infinite words, 
called the shift-space of X. 

Give X the discrete topology. Thus E inherits the product topology (i.e. the topology 
whose basis is sets of the form YViLi such that there is some M with Ai = X for i > M). 
We define the maps (7^ : E — >■ E by ai{w) = iw, for all w g E. We also define the map 
(7 : E -> E by a{wiW2'W3 ■ ■ ■) = W2W3 . . ., for wiW2'W3 . . . € E. 

Definition 2.1. Let Khea compact metrizable space, let X be a finite set for i*"; : K ^ K 
continuous injections such that K = Ui(^xFi{K). We call the triplet C = {K, X, {^iligx) 
a self-similar structure on K if there is a continuous surjection tt : E — > /-C such that 
Fi o TT — TT o ai, holds for all w eX* . 

Ifw — W1W2W3 . . . Wr, then we define Eui :~ F^-^ o F^^ o • • • o F^j^ and = F^ (K). 

Definition 2.2. Let C = {K, X, {^i}jgx) ^ self-similar structure on K. The critical set 
of £ is defined 

ciC)^TT-n y {K,nK,) 
and the post-critical set 

V{C) = U a"{C) , 

n>l 

where cr is the shift operator on X~'^. We abbreviate the critical and post-critical sets by C 
and V when there is no danger of confusion. A self-similar structure is called post critically 
finite (p.c.f.) if T-" is a finite set. 

We shall define the fractal boundary to be Vb := 7r(7'), and recursively define Vn — 
Uu,GX"-F«i(Vo), in particular /\„ Ci Ky (- Vn for any u ^ v E X". 

For a finite set V, define the vector space i{V) := {u : V M.} with pointwise addition 
and scalar multiplication. i{V) is the set of vectors indexed by V. 

We use equation 1 1.11 to build a Dirichlet form on K from forms on Vq. We define the 
conductance matrix C whose rows and columns are indexed by Vq, to be non-negative 
symmetric with diagonal entries which are identically 0. That is Cpq > 0, Cpq — Cqp, and 
Cpp = for all p and g S Vb- We define a bilinear form 4) : ^(Vb) x £{Vo) M by 

^o(u,w) = ^ ("(P) (i-b) - (9)) Cp, 

For and bilinear form, we define the induced quadratic form (^o{u) := <S'o{u, u). 

For p = {po, . . . , pn-i), inductively define an energy form on T4, by a map 'I'p which 
takes an energy from onVn-i and produces an energy form on Vn by 

N 

Sn{u,v) p{Sn-i){u,v) := '^PiSn-liuo Fi,V O Fi). 
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We define the energy form S" on K 

S{u,v) = lim (^„(m|v„,w|v„). 

and refer to set of u such that S{u) is well defined as the domain of 

However, we can not chose p and C arbitrarily, we require some consistency between 
approximating energies Recall Vn-i C Vn- Define the operator Tr which takes a 
Diriclet form on Vn and produces an energy form on Vn-i- For u E £{Vn-i), 

TrfS'„(u,w) inf {S'nig) \ 9 e £{Vn), g|v„_i = u} 

The consistency condition we need is that TtS'i = S^. If this is the case, we can find an 
element u e (-{Vn) with {t|v'„_j — u and (?„({t, w) — Tr(?„_i(u,M) — (?„_i(u,m). u is 
called a harmonic extension of u, we shall discuss such extensions more in section[6l 

Determining a harmonic structure on a self-similar structure K is equivalent to finding 
a connected matrix C (We call an n x n matrix M connected if for alH, j — 1, 2, . . . , n, 
i 7^ j, there is i = ii, 12, . . . , u- = j such that Mi^i^^^^ 7^ for £ = 1, 2, . . . , fc — 1.) and p 
such that Tro^'p((?o) = <§o- 

3. Energy Forms as Matrices 

For no < ni, we think of R"" ® R"° as tiq x no matrices, M"i (g) R"i as n-i x ni 
matrices. Consider the collection of ni x np -matrices {V^i}^]^ with entries which are all 
zero except for exactly one 1 in each row and at most one 1 in each column. Thus {^/'i }^]^ 
are linear injections V'i : ^ R"^ . Assume that 

spanlV-^lR"") I i = l,2,...,iV} = R"i 

Define the map 

* : (0, 00)^ X R"« «) R"o ^ R"i (g) R"i 

by 

N 

= ^p,V^M^• 
We can also think of ^I' as defined on simple tensors by 

N 
i=l 

and extended linearly. 

Elements of R"^ (g R"^ can be written as block matrices 

,r f A B 
^^=U D 

where A e R"" (g) R"o, B e R"" (g R"i, C e R"i (g) R"^^ and D e R"i-«o (g R"i-"o. 
We define the map 

Tr : R"i g) R"i ^ R"» ® R"« 

by 

Tr(M) = A - BD^^C. 
This is called the Schur complement in general. Finally, we define the map 

A = Tro^- : (0, 00)^ x R"° (g R"° ^ R"° (g R"". 

In the fractal theory, Vb ) is R"« and ^ ( Vi ) is R"i , where no = | | and ni = | Fi | . 
is the map which extends a function in £{Fi (Vq)) to a function ^(V^i) with values identically 
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on \ Fi{Vo). There is a natural correspondence between bilinear forms and matrices 
given by mapping a matrix M to a form 

S'Mix,y) = {Mx)-v = y'^Mx. 

The energy forms in the previous section correspond to matrices M that have non-negative 
off-diagonal entries and satisfy Ml — where 1 is a vector with all entries equal to 1. In 
particular, Cpq are the off diagonal entries. 

It is simple to see that — S'-qi^p j^.i^ . It is shown in | KigOl) chapter 2 that Tr <§m = 

<oTi M, where Tr <Sm is defined as in the last section. 

Thus finding a self-similar energy form is the same as finding p e (0, oo)^ and a con- 
nected matrix M e K"" ® E"" , as above, such that 

(3.1) A{p,M)^M7 

4. Some basics of electric circuits 

The techniques used in computing the harmonic structure on the fractalina gasket come 
from the basic theory of electrical networks. For a full introduction to these networks and 
their relationship with Dirichlet forms see liDS84l . 

Consider V, a finite collection of nodes and a symmetric matrix of conductances C = 
{Cpq}^ qey ^i'^h positive entries. We often draw F as a graph, where there is an edge 
between p and q if Cpq > 0. Likewise if p is not connected to q in this graph, then Cpq — 0. 

Think of Cpq has being the strength of an electrical connection between p and q. On the 
other hand, one can consider the resistance between two points: Rpq = l/C'pq if Cpq > 
and Rpq = oo if Cpq — 0. If we take a function u : y ~> M then we can take the energy of 
the state 

p.qev 

We are most often interested in computing the effective resistance between points p and q in 
V. Informally, this is the resistance if one node was hooked to the positive end of a battery, 
while the other was hooked to the negative. Formally, is defined by 

„ ] , = inf {^(u, u)\u:V^R, f(u) - 1, f(u) - 0} , 

fox p ^ q and Rctf{p,p) — 0. It can be shown that i?eff is a metric on V (see |Kig01| ). 

When is a self-similar resistance form, eff'ective resistance scales with the self-similarity 
of the fractal. Thus if we can calculate the effective resistance between two points in the 
fractal on two levels of levels approximations, then we can obtain the resistance scaling 
factor. This is the technique we employ to calculate the scaling factor of the fractalina, and 
is the same strategy used in llBCF+0 71 to calculate the resistance of fractal A^-gaskets. 

This reduces the problem of finding the resistance scaling factor to that of finding the re- 
sistance of an electrical network. We shall need two tools from the basic theory of electrical 
circuits to reduce the conductance networks — Kirchhoff's laws and the A-Y transform. 

Kirchoff"s two laws state that connections in parallel add conductances and connections 
in sequence add resistance. This allows us to reduce the graph as follows, if there are two 
connections between nodes p and q with conductances Cpq and Cpq, then the total conduc- 
tance is Cpq + Cpq . Likewise if there is a connection between node p and r with resistance 
Rpr and between nodes r and q with resistance Rrq, and there are no other connections to 
r (i.e. Crs = for s e V with s ^ p and s ^ q), then we can combine the two connections 
with a connection between p and q with resistance Rpq — Rpr + Rrq- 
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Figure i . The A-Y transform. 



If we can reduce a graph to one connection between p and q, then the resistance of this 
connection is the effective resistance between these points. In particular, if we reduce the 
graph to be such that there is only one path between p and q which does not cross itself, 
then the effective resistance between p and q is the sum of the resistances of each edge of 
this path. 

The other main tool we use in reducing networks is the A-Y transformation. This tells 
us the process by which to transform a triangle network to a network which looks like a 
"Y" in a way that the effective resistance between any two nodes remains the same. If the 
resistances are set up as they are in figure [T] then 

RabRac , r> + ^&''c + ^a^c 

and Rab = ■ 



Rab + Rbc + Rac 

This transformation is useful because it is easier to find the effective resistance between two 
nodes of a Y network than a A network. For example, i?cff (a, b) ^ Va + ri,. 

5. Fractalina Fractal 

Fractalina was introduced by Robert Devaney using standard techniques in iterated func- 
tion systems IICDF98II . The construction of the fractalina fractal is similar to that of the 
Sierpinski triangle. 

Let Pi, P2, and P3 be vertices of an equilateral triangle in the complex plain, fractalina 
is the unique compact Hausdorff set invariant under the iterated function system (IPS) con- 
sisting of 

Fi{z) = F2{z) = and F3 = . 



Similar to the IFS of the Sierpinski gasket the contractions Fi have contraction factor of 
1/2, however, the contraction F3 includes a rotation around P3 of 180 degrees. 

If assumptions are weaken so that Pi are the vertices of any triangle, then the maps above 
result in a homeomorphic fractal. We assume equilateral, to allow us to make deductions 
based on the symmetry of the resulting fractal. 

We assign resistances between each vertex of the Level 1 fractal, the graph hour-glass 
graph on the left in figure|2l In addition to the fractal boundary, the set {Pi, P2, {Pi ) , F3 {P2 ) }, 
we include P3 in our first level network, because removing this point disconnects the fractal. 
That is we assume 

Vo = {Pl,P2,P3,F3iPi),F3{P2)}. 
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Figure 2. Level 1 and Level 2 of fractalina 




In addition, we assume the network consists of two triangular networks, {Pi, ^2,^3} and 
{F3{Pi) , F3{P2) , P3} ■ The second level is attained by applying the contractions, and we 
end up with the graph on the right in figure |2l 

Theorem 5.1. The resistance scaling factor between any two levels of the fractalina is 
(3 + ViT)/16. 

Let k represent the scaling factor associated with the Fractalina contractions. Since 

F3({Fi, P2, P3}) = {F3{Pi),F3[P2), F3} , 

the resistances between elements on the right hand side are are those of {Pi , P2 , ^3 } scaled 
by k. Noting the reflection symmetry through the vertical axis, we arrive at the configu- 
ration of resistances found in figure [3] This figure also shows the application of a A-Y 
transformation to the two triangles. The three new resistances in figure[3]given in terms of 
i?i, R2, and k by 

fcPlP2 « k{R2f + {R2f ^ R1R2 

a ~ , n — , and 7 = . 

2P2 + P1 2P2 + P1 2P2 + P1 

We now consider the level 2 fractalina. Figure |4] shows the unsealed resistances between 
vertices. 5 applications of a delta-Y transformation yields resistances according to figureHl 
a' = 6 — 'y,S = ri = f3,e = 2a, and l = 2j. The resistance C is equal to a, but this is 
unimportant as we shall not be considering the eff'ective resistance for any pair of points that 
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a path between which would require these connections. A restructuring of the resistances 
using Kirchhoff's laws gives the lower left side of figure |4] with the new resistance 

^ 2kRiR2 + kRl + {R2Y 

K — la ~\- p ~ . 

2i?2 + i?i 

One more A-Y transform results in the resistances 

(2fci?li?2 + fc(i?2)' + {R2?? + + (i?2)')(4fciT'iiT'2 + iKR^f + 2(^2)' + IR1R1) 

(2i?2 + Rl){ikRiR2 + 2/c(i?2)2 + 2i?2 + 2R1R2) 

{2kRiR2 + fc(i?2)' + {R2?){2R2Ri) + {RiR2){ikRiR2 + 2fc(fl2)' + 2(i?2)' + 2flifl2) 

(2i?2 + i?l)(4fci?ii?2 + 2fc(i?2)2 + 2(i?2)2 + 2R1R2) 

By comparing the effective resistance between Pi and P2 on both levels, we get that 
k = 27/27'. Comparing the resistances between (Pi ) and Pi , we get that k = (a + /? + 
7)/(q;' + + 7'). So now there are two equations and three unknowns: k, Ri and i?2- 
We can reduce the number of unknowns by assuming, without loss of generality, R2 = \. 
We can do this because resistance is a quadratic form, and we are only interested in the 
relationship between Ri and R2 up to a constant multiple. This is manifest in the fact that 
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all the above resistances are homogeneous rational functions in i?i and i?2- We reduce the 
first equation 



^ 2kRi + k + l + R, 
^ ' AkRi +2k + Ri 

Solve for i?i in terms of k, 

1 - fc - 2fc2 



Ri 



4fc2 - fc - 1 



Assume fc = (3 + v4T)/16. This implies that Ri = (-1 + \/4T)/4. Of course, it must 
be verified that these values satisfy the second equation, which it also satisfies. Thus the 

scaling factor of the resistance form on the fractalina fractal is ^-i . 

^ 16 



6. Pillow Fractal 

The pillow fractal is a self-affine fractal in that the contraction mappings that produce 
this fractal are not rigid. The basic building blocks of the pillow are rectangles and there are 
two contraction maps. Both contraction maps rotate the fractal 90 degrees, gluing together 
two of the corners while leaving the edges separate. FigureH] shows this construction. 

In 0NTO81 . the pillow is constructed as the limit space of a self-similar group. A recur- 
rence relation is determined from the action of this group, which gives equations to deter- 
mine the resistance of the Schreier graphs approximating the limit space. These Schreier 
grpahs are different than the graphs used used in the current paper to approximate the pillow, 
however, the results are consistent. 

The energy on the pillow is built up from approximating graphs, the 0th level is a "rec- 
tangle" in that it has four nodes Fq — {vi,V2, ws, W4}, a node at every corner. We are 
assuming rectangular symmetry of the pillow. In particular the energy form is invariant 
under horizontal and vertical flips as well as 180 degree rotation. More precisely, for the 
maps cr = {viV2){v3V4) , T — (wiW3)(w2t'4), P — (wit'4)(t'2i'3) in the symmetric group of 
Fo, 

^o(/) = S'oif o cr) = S'o{f o r) = S'oif o p) 

This invariance implies that we can assume that C^-^v^ — Cy^v^ ■= Ci, C'unj2 = C'i>3U4 '■= 
C2, and C^i^4 — C^2f3 ^3- 

Define Fi = _Fo U {ui, U2}, and the first contraction maps vi — ?> vi, V2 — > U2, V3 -> V2, 
and 114 — >■ U2, while the second maps t^i — > vi, V2 — >■ U2, V3 — > V2 and V4 — ?> U2. The 
conductances of Fi are determined by these contractions, noting that we get two edges 
connecting ui to U2 with conductance Ci (equivalently, we can think of C,jj,j2 — 2Ci). 
Figure |6] shows Fq and Fi. 

The rectangular symmetry from Fq can be extended to Fi, that is that S2 is invariant 
under (7 = {viV2){v3Vi){uiU2), f = {viV3){v2V4), p = {viV4){v2V3){uiU2). Notice that, 
for example alp^ — a 

Theorem 6.1. The resistance scaling factor between any two levels, Si and <^i+i is 

The rest of this section is the proof of this. To find the resistance scaling for the pillow, 
we employ a different strategy than in section |5] We look at specific states of the fractal, 
and see how they correspond to the graph at each level by assigning potentials to each node. 
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Fo 

I 7\ 



1 i 1 

Fi 

Figure 5. Contrustion of Fi from Fq 



To compare states consistently between graphs, we take a function / : Fq ^ K, we get 
a function on Fi by harmonic extension — the function / such that f\pg — f, and 

(6.1) E ^ 

for all verteces p ^ Fq, and Cpq is the conductance between p and q. The harmonic exten- 
sion / is the energy minimizing extension of /. That is to say S'i{f) < A (5) for any state 
g such that glp^ = /, or that A (/) = Tr A (/). For more on this see IIStr06llDS84t . 

We make some important observations about harmonic extension that are straight for- 
ward to verify. First, harmonic extension is linear: if / and g are the harmonic extensions 
of / and g respectively, then f + eg if the harmonic extension of / + eg, for any constant 
c. Second, the extension of a constant function is constant. This is easy to see, because 
the energy of a constant function is 0, so the constant function is an energy minimizing 
extension. Finally, harmonic extension is invariant under symmetries of Fq and Fi, that is 
to say if / o (7 = /, then f o a — f (similarly for r and p), this can be seen by examining 
the definition in l6.1l 

Comparing the energy of the state / to the energy of / will give the resistance scaling 
factor, S'o{f)/S'i{f). However, this calculation leaves unknown quantities. We will have 
to consider three functions /i , /2 , /s : Fq — > M to solve for these quantities, and obtain the 
scaling factor. 

Configuration 1 Define fi{vi) = fi{v2) — fiivs) — and /i(w4) = 1 
To ensure that we get the appropriate energy we calculate the energy minimizing poten- 
tials at the unaccounted for verteces of Fi, = x and /i(u2) = y as in figure[2l The 
values X and y satisfy equation [6.1l Written in terms of Ci, for i = 1,2,3, 

2Ciy + C2 . 2C\x + C3 

(6.2) X = —— and y — 



2(Ci + C72 + C3) ^ 2(Ci + C2 + C3) ■ 
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Figure 7. Configuration 1 



These equations are used to calculate the energy form for the first and second orders, whose 
ratio is 

^o{fi) ^ Cl + C2 + C3 

^i(/i) C-iCl + 2(a; - yf) + C2{x^ + 22/2 + {x - 1)2) + C3(2a;2 +y^ + {y- 1)2) " 
We complete similar calculations for two more cases. 

Configuration 2: Define /2(v2) — /2(f3) = and /2(t'i) = /2(f4) = 1- Notice that 
/2 o P = /2, so /2(w2) = /2(/5(wi)) = /2(mi). Furthermore, /2 + /2 o r is the constant 
function 1, so /2(wi) + /2(f («,i)) = 2/2(«,i) = 1. Thus f{ui) = 72(^2) = 1/2. 
Since we know all the values of /2, we calculate 

Mh) ^ 2Ci + 2C2 
A(/2) 2C1 + C2 + C3 

Configuration 3: Define fsivi) = 73(1^2) = and faivs) = faiv^) = 1. Similar to the 
computation of /2, /3 o cj = and /3 + /s o r = 1, so /3(ui) = /3(u2) = 1/2. Leaving 
us to compute 

Mh) 2C1 + 2C3 



^l(/3 



C2 + C3 
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s 

Assume — — 2^1'^ for the latter two equations, and without loss of generality, let C\ = 
1. This implies that 

2l/3 ^ 2 + 2I/3 = ^ + ^^3 

2 + C2 + C3 C2 + C3 

Solve for the two conductances we get that C3 = 2^^/^ and Ci = 2^^/^ + 2^/^. It now 
must be verified that these values satisfy the equation for (§b(/i)/<ft(/i)- Apply the values 
of Ci to the equations 16.21 we get that x = 2^^/'^ and y = — 2"^/'^. Applying these 

^0 ( fi ) — 1 / 

results into the equation above we get that ^ = 2 ' which proves consistency. 
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